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Abstract 

We average the decay width of a heavy meson in the t' Hooft model over 

the heavy quark mass M with a smooth weight function. The averaging has 

support over a few resonance spacings. We use the previously determined 

heavy meson decay width which differs from the free quark width by a 1/M 

correction. In contrast, we find that the averaged meson and quark widths 

differ by a 1/M 2 correction. We speculate on the relevance of our results to 

the phenomenologically relevant case of 3 + 1 dimensional QCD. 
PACS number(s): ll.10.Kk, 11.15.Pg, 13.25.-k 
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Quark-hadron duality has been part of the lore of strong interactions for three decades. 
Bloom and Gilman [3],Q] (BG) discovered duality in electron-proton inelastic scattering. 
There, the cross section is given in terms of two Lorentz invariant form factors W\ and 
W2 which are functions of the invariant mass of the virtual photon, q 2 , and the energy 
transfer to the electron, v. Considering the form factors as functions of the scaling variable 
uj = q 2 /2Mu, they compared the scaling regime of large q 2 (and large v) with the region of 
fixed, low q 2 . They determined that, for each form factor, the low q 2 curves oscillate about 
the scaling curve, that identifiable nucleon resonances are responsible for these oscillations 
and that the amplitude of a resonant oscillation relative to the scaling curve is independent 
of q 2 . Moreover, they introduced sum rules whereby integrals of the form factors at low and 
large q 2 agree and noticed that the agreement was quite good even when the integration 
involved only a region that spans a few resonances. 

Poggio, Quinn and Weinberg || (PQW) applied these ideas to electron-positron anni- 
hilation. While BG compared experimental curves among themselves, PQW compared the 
experimental cross section to a scaling curve calculated in QCD. They noticed that the 
weighted average of the cross section a(s), 



is given in terms of the vacuum polarization of the electromagnetic current with complex 
argument, 



and argued that one can safely use perturbation theory to compute this provided A is large 
enough. This procedure was better understood with the advent of Wilson's [|J Operator 
Product Expansion (OPE). It is interesting to point out that the prediction of PQW based 
on the two generations of quarks and leptons known at the time did not successfully match 
the experimental results. When PQW allowed for additional matter they found a best 
match if they supplemented the model with a heavy lepton and a charge 1/3 heavy quark, 
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FIG. 1. Normalized e + e~ cross section averaged over the squared center of mass energy s as 
in Eq. (Q) with A = 3GeV 2 . The dashed (red) curve is from PQW ||. The solid (green) curve is 
ours and shows the effect of including the narrow charm-onium resonances. 

anticipating the discovery of the tau-lepton and b-quark. It is also interesting that PQW 
did not include in their average of data the contribution of charm-onium resonances. When 
this is done the average cross section is raised significantly at low s, leaving the higher s 
region unaffected, as shown in Fig. [p. 

In an attempt to understand the origin of quark-hadron duality we have computed both 
the actual rate and its "scaling limit" from first principles in special situations. In Ref. || we 
computed the semi-leptonic decay rate and spectrum for a heavy hadron in the small velocity 
(SV) limit. We showed that two channels, B — > Dev and B — > D*eu, give the decay rate 
to first two orders in an expansion in 1 /m& and that to that order the result is identical to 
the inclusive rate obtained using a heavy quark OPE as introduced in Ref. |J. The equality 
holds for the double differential decay rate if it is averaged over a large enough interval of 
hadronic energies. The computation demonstrates explicitly quark-hadron duality in semi- 
leptonic 5-meson decays in the SV limit, but really sheds no light into the mechanism for 
duality. In particular, it is puzzling that duality holds even if the rate is dominated by only 
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two channels. 

More recently we attempted to verify duality in hadronic heavy meson decays. In Ref. |7j 
we considered the width of a heavy meson in a soluble model that in many ways mimics 
the dynamics of QCD, namely an SU(N C ) gauge theory in 1 + 1 dimensions in the large N c 
limit. This model, first studied by 't Hooft exhibits a rich spectrum with an infinite 
tower of narrow resonances for each internal quantum number, making the study of duality 
viable. We considered a '5-meson' with a heavy quark Q and a light (anti-)quark q of 
masses Mq and m, respectively, which decays via a weak interaction into light qq mesons. 
To leading order in 1/N C the decay rate is dominated by two body final states: if TTj denote 
the tower of gg-mesons, the total width is given by T(B) = — > Tij-Kk), where the 
sum extends over all pairing of mesons such that the sum of their masses does not exceed 
the B mass, fij + fik < Mb- The main result of that investigation was that there is rough 
agreement between T(B) and the decay rate of a free heavy quark, T(Q). When considered as 
functions of Mq the quark rate is smooth but the meson rate exhibits sharp peaks whenever 
a threshold for production of a light pair opens up. This is due to the peculiar behavior 
of phase space in 1 + 1 dimensions, which is inversely proportional to the momentum of 
the final state mesons. Nevertheless, in between such peaks it was found that the relation 
T(B) = r(Q)(l + 0.14/ Mq), in units of g 2 N c /n = 1, holds fairly accurately. 

In this brief paper we consider the effect of local averaging on the results of Ref. 0. The 
main result is that when averaged locally over the heavy mass Mq the agreement between 
T(B) and T(Q) is parametrically improved. In fact, for the averaged widths we find 
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Remarkably, the correction of order 1/Mq has disappeared. 

The computation uses the numerical results of Ref. J7|. The averaging is defined by 

r(M = J -&% — . An2/ 2 V ; . 4 

J ^min 

The limits of integration are the lowest and highest heavy masses available from Ref. J7|. 
The width a was taken to be a = 1, the scale of the strong interactions in our units. The 
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FIG. 2. Log-difference of averaged meson and quark decay rates as a function of log(M(g) 
(dashed lines), and straight line fit (solid line). The different dashed lines correspond to powers 
n = —2, —1, 0, and 1 in the average of Eq. (|j) (n = —2 is highest at small Mq). The straight line 
is -1.71og(M Q ) + 1.45. 

integer n was varied between —2 and 1 to study the effect of emphasizing low or high masses 
in the average. The result is shown in Fig. |2| where log[(r(5)) — (T(Q))] is plotted versus 
log(MQ). The peculiar behavior at low and large Mq (outside 1.8 < logMg < 2.3) is due 
to the finite endpoints in the integral defining the average as can be readily checked by 
averaging simple smooth curves. The nearly straight intermediate region is therefore what 
concerns us. The straight line in Fig. ^| is a linear fit by eye, which is accurate enough since 
there is some variation in the results depending on the value of the power n used and what 
range of log Mq is fitted. We have checked that small variations in the linear fit do not alter 
our conclusions. We then find the best fit of the exponential of this line to the function 
a + b/M Q + c/Mq, yielding a = —0.01, b = 0.47 and c = 4.69. The value of a is consistent 
with zero, and the quality of the fit is improved by dropping this degree of freedom (and 
gives Eq. <g)). 



Since, when plotted as functions of Mq, T(B) is (almost) always above T(Q) and the 
difference seems constant except at the narrow peaks of T(B), the question immediately 
arises as to how the averaging procedure can turn the constant difference into one that 
decreases as 1/Mq. The answer is suggested by the plot in Fig. [I]. Recall the average curve 
of PQW does not include the effect of the narrow charm-onium resonances. Our curve 
includes these narrow peaks in the averaging and when compared to the curve of PQW it 
is enhanced at low s. The same is apparently occurring in the averaging of T(B). The area 
under the narrow peaks at low Mq is larger than at high Mq, tilting the curve slightly. 



The result of Ref. Q was criticized in Ref. |§ (see also Refs. ||10| - |12|). The inconsistency 
between the analytic results of Ref. and the numerical results of Ref. |7j was blamed on 
numerical inaccuracies of the latter. However, it is clear now that this is unlikely. It seems 
impossible to understand what error, random or systematic, would conspire to change the 
behavior of the difference between widths from constant for non-averaged widths to decaying 
as 1/Mq for averaged widths. 

The main result of Ref. JJ, that there are no 1/Mq corrections to the partonic width, is 
easy to derive. Consider the decay Q — > q + q'q' where the weak interaction giving rise to the 
decay is from the operator cpi^Q q'^^q' , and take the mass of the q' quarks to vanish. In this 
limit the current q'^f^q' couples only to the lowest meson in the tower of q'q' states, which 
is massless. This follows from current conservation: from (0\q' , y^'~f5q'\p) = fp^ it follows 
that fp 2 = 0, so / 7^ only if p 2 = 0. It follows that the heavy meson decays into pairs 
of the form 07r n where <fi is the massless q'q' meson and 7r n are the mesons in the tower of 
qq" states, where q" is the (light) spectator quark. The amplitude for B — > <j)7r n is therefore 
determined by the B — > 7i n form-factor at momentum transfer q 2 = 0. The calculation of 
T(B) is therefore identical to the semi-leptonic width at zero ev invariant mass. This is 
guaranteed to have no 1/Mq corrections by the argument of Ref. ||. 

Ref. attempts to extend this result to the case m! ^ 0. This is a very delicate 
matter and we disagree with the procedure presented there. In particular, the perturbative 
calculations are out of control. The dimensionless expansion parameter in the 't Hooft model 
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is g 2 N c /7cm 2 , where m is the smallest quark mass in the problem. Take for example the 
vertex correction to the qQ current, and let k be the momentum carried by the current. If 
k 2 = identically then, as shown in Ref. ||, the vertex correction vanishes. However, for 
k 2 7^ the leading term in the vertex correction at one loop is proportional to g 2 N c /i:m 2 
(see Ref. [0 and, in particular, appendix B of Ref. [0]). Moreover, in the limit k 2 — > 
the vertex correction is non-vanishing. In fact, if all light quarks have mass m, the leading 
contribution to the heavy quark width at one loop is (after a lengthy computation) 

r(g)( i-ioop ) = 3 ^^^ r(Q)(trec) _ (5) 

It is dangerous to study the case of non- vanishing light quark masses as perturbations about 
vanishing mass. 

The question immediately arises as to why is T(Q) a good approximation to T(B) in 
the first place, since perturbation theory clearly breaks down at small m. Operationally 
the answer to this question is that the re-summation of the gluon ladder that gives the 
quark form factor of a vector current is well approximated by the free quark form factor. So 
although perturbation theory looks hopeless, the re-summed vertex is almost like free. 

We speculate on the relevance of our results to the phenomenology relevant case of 3 + 1 
dimensional QCD. The heavy meson width is not expected to have singular peaks (as a 
function of Mq) as in the 't Hooft model. However, it is entirely possible that it oscillates 
about the quark-dual rate. We cannot rule out, and indeed this work suggests it is entirely 
possible, that the amplitude of these oscillations decreases with Mq only as one inverse 
power. Once averaged over Mq with a weight function of width of order 1 GeV the rate 
may display oscillations that decrease faster with Mq , say, as 1 / Mq , as apparently indicated 
by "practical OPE" arguments |J. However, it is impractical to average over heavy quark 
masses, thus leaving us with uncertainties of order 1/Mq in our predictions of hadronic 
widths. Some evidence for this was presented in Ref. |15|] where it was observed that the 
6-quark width agrees better with experimental hadronic widths if the quark mass is replaced 



by the B or masses, respectively. In a similar vein, Ref. [0 shows how 1/M violations 



to local, but not global, duality may occur in £?-meson correlations. 

In summary, we have shown that quark-hadron duality in heavy meson decays in the 
't Hooft model is accurate only to order 1/Mq, but the accuracy is promoted to order 1/Mq 
for local averages over Mq of the widths. We speculate that the same effect occurs in four 
dimensional QCD. Since averaging over 6-quark masses is impossible, it is safe to assume 
that the heavy meson and baryon widths are at best computed with error of order 1/Mq. 
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